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Abstract 

We account for the interaction of the Bose-condensed fraction with the 
normal phase in an effective dynamical equation such as the Gross-Pitaevskii 
equation. We show that the low-energy excitations can be treated as sound 
waves with speed dependent on the condensate density. This allows us to re- 
duce the problem to the calculation of the determinant of the Laplace operator 
on a curved space and apply standard methods of quantum gravity to get the 
leading logarithmic contribution of the determinant. This produces the first 
quantum correction due to the noncondensed fraction to the Gross-Pitaevskii 
equation for the condensate. The correction describes an additional quantum 
pressure in the condensate and evaporation-condensation effects. 
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1. Last years there was a new growth of the interest in Bose-Einstein condensa- 
tion (in particular, for spatially inhomogeneous situations) due to the recent success 
of the experimental observation of Bose-Einstein condensation for systems of spin 
polarized magnetically trapped atomic gases at ultra-low temperatures [l], |], |3| and 
further investigation of their collective properties ||. These experimental results 
stimulated the development of a theory of a nonuniform Bose condensate and its 
collective excitations. Moreover, new exciting problems such as a description of the 
evolution of the Bose condensate from a relaxed trap, stability and dynamics of 
the collapse of the condensate for Li 7 atoms, collective excitations, heating-cooling 
phenomena, various coherence effects for the condensate and so on are being posed 
both theoretically and experimentally. 

All the problems listed above require an accurate account of the interaction 
of the condensate with the non-condensed fraction. Formally, the problem can 
be tackled using the Bogoliubov-Popov equations which represent a set of coupled 
linear and nonlinear Schrodinger equations; the nonlinear Schrodinger equation (for 
the condensate fraction) coupled via a potential term with a pair of Schrodinger 



equations (for the excited modes) whose potentials, on their own, depend on the 
solution of the nonlinear Schrodinger equation. To find an effective equation for 
the condensate we need to resolve a pair of Schrodinger equations with arbitrary 
potentials and find the solutions of the equations as functionals of the potentials. It 
is obvious that this is hardly a realistic task to perform. 

However, as we show in the paper, it is possible to calculate the corrections 
from the non-condensed fraction and, surprisingly, they look very simple in specially 
chosen notations. The answer comes from exploiting a technique originally developed 
to solve problems in quantum gravity, which seems to be quite far from the Bose 
condensation. To make use of it, we note that the low-energy modes of the non- 
condensed fraction can be considered as sound waves whose velocity depends on 
the density and velocity of the condensate. This shows that the interaction of the 
modes with the condensate has the same nature as the interaction of matter fields 
with the background gravity field. In this framework the Gross-Pitaevskii equation 
is an analogue of the Einstein equation where the energy-momentum tensor is due 
to presence of quantum matter fields, that is excited modes. It is not surprising 
then that to calculate the correction we use methods from quantum gravity. More 
precisely, we use functional integration, (- function renormalization for determinants 
of elliptic operators and the Seeley expansion of the heat-kernel for the Laplace 
operator on a curved space in an external field. The final form of the correction to 
the free energy of the condensate, written in terms of the curvature tensor then has 
a very simple and manageble form (]g). 

2. We are looking for the effective action which describes all physical quantities 
in the system (for an easy introduction to the formalism of the effective action see, 
for example, Appendix A of Ref. 0). For example, the effective action provides all 
of the Green's functions in the same approximation used to calculate the effective 
action itself. The equations of motion for the average values of field operators are 
given by the equations for a minimum of the effective action. We obtain them in a 
hydrodynamic one-loop approximation. To clarify the description the hydrodynam- 
ical variables, density and velocity, should be used instead of microscopic variables. 
More precisely, we start with the action 

+ (v- + K4> + 4>f 

and change field variables to the hydrodynamic ones: 

ip(x, r) = y/p(x,T)e-^ , il> + (x, r) = y/p(x,T)e^ . 
This transforms the action ([[]) into the form (up to a complete derivative term): 
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It is not difficult to see that the classical equations of motion for the action lead to 
the equations: 

dc ( c 2 2 \ 

^ + yy- + 2v-2 l i + Alp-— V 2 v/pJ =0, (3) 
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| + V(cp)=0 
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which are equivalent to the Gross-Pitaevskii equation || after the introduction of the 
velocity variable c = — 2V<p instead of p. These equations look like hydrodynamic 
equations for an irrotational compressible fluid which is the condensate superfluid. 

Now let us shift our variables in ([I]) by the zero-order mean field solution (the 
condensate amplitude) to find the one loop quantum correction for the effective 
action: 

ip -> ip + X , ip + + x + 

and keep only quadratic terms of \ an d X + since it can be shown || that this is 
the only relevant contribution for the case of trapped gases. Then the action (|TJ) 
transforms to S(ip, + Si(ip, X> X + ) where the term Si can be written as: 



Si(ip,ip + ,X,X~ 



dr jdvL X + ^ 



Vx + Vx + {p- (V^) 2 + 2lp} X + X 



-le 2i ^ 



PX X 



le~ 2i * 
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In the last equation the potential v was excluded using the mean field equation of 
motion neglecting the kinetic term. The Hamiltonian corresponding to S\ is 



H= dV 



Vx + Vx + {p~ (V<p) 2 + 2lp}x + X - le 2i *px + X + - le~ 2i<p pxx 



Now we perform the following canonical transformation: 



X = y/fxf 
The Hamiltonian takes the form 



a 

— + ia 
2p 



X + = VP? 
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H = JdV 



p(Va) 2 + ppa 2 + {l + f-}a 2 + a{VpVa) - d(V^Va) + \Apaa 



4p 



(Vp) 



-« 2 + ^(Va) 2 + 



(Vp) 2 ^ 2 
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Last eight terms vanish in the hydrodynamical limit since p is a large variable 
(V<p ~ y/Tp and (p ~ Ip), so we have 



H= dV 



p{Va) 2 + ppa 2 + {l + j^}a 2 + a(V<pVd) - d(Vy?Va) 
The corresponding action (in the Feynman functional integral) is 



Si 



^drfdVp 
-oo J or 



p{Vaf + ppa 2 + {l + -^}a 2 + 2a(VpVa) 



Integrating the functional integral over the field a we get the action for a field only: 



Si = /"dr JdV j-p(Va) 2 + /(^ -2VpVc 



a I — ppot 



where / = ^p~n ■ The corresponding determinant from the integration does not 
contribute to the effective action at this order of the perturbation theory. 

Now we will once again use the existence of the large parameter p. Let us 
introduce the following large number p = max{p} such that p — ~ 1, and new 

space-time variables t = 4Zp r and y. t = ^4lp Xi. The action Si can be represented 

as 



Si 



1 



A^d^ad^a + pipe? 



ptpa 



dt dV 



where 



v = 2V<p , / = Alf 



and the matrix A has the form 
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( -f fVl fj&2 f?3 \ 

fyi p-fvi -fviy 2 -fyiv 3 

fp2 -fpiV 2 P-fvl -fv 2 V3 

V /«3 ~fVlV 3 ~fv 2 V 3 p-fvi J 



Our next step is to cast the action in the covariant form. To do this we introduce 
auxiliary metric g^ u such that 



-fiu 



-det (II AH 



/-det(II^H) 
This leads to the form for the metric Hg^H: 

/ -p~/f + v 2 Vi v 2 v 3 \ 



\9nu\ 



\ 



Vl 

v 3 



In this metric the action has a covariant form 

1 



where 



E 



Si 
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dy 
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aAa + Eot 
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Now the effective action is Ti = — |tr ln[(64po/ 3 ) 1 ^ 2 (A. + E)] with Laplace operator 

At this point the quantum gravity analogy comes into the game and let us 
evaluate Ti using ^-function regularization |7|] for |tr hi[(64poZ 3 ) _1 / 2 (A + E)\ after 
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performing the Wick rotation to make the operator elliptic. In this regularization 
the quantity is finite and has the following property: 

tr ln[(64p / 3 r 1/2 (A + E)} = tr ln[A + E} — -tr ln(64p / 3 ) ($o(A + E) - L(A + E) 

where $o(A + E) = Jdydr ^—g^/ (A + E) is so-called zeroth Seeley coefficient 
and L is a number of zero modes. For our goal is important to realize that in the 
regularization tr In [A + E] is of order of the zeroth Seeley coefficient which, on its 
own, is of order unity. 

The existence of the large parameter in the system allows us to make the next 
step. Indeed, the multiplier (64p ^ 3 )^^ 2 gives us the possibility to express the main 
contributions to the determinant such as 



tr ln[(64poZ 3 r 1/2 (A + E)] 



itr ln(64p / 3 ) 



dxdt(^ (A + E) - L(A + E) 



The Seeley coefficient \l/o f° r the Laplace operator in 4D is well-known and can be 
found, for example, in Ref. 0: 



(47T) 
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where R, R pi/ and R^ uap are the scalar curvature, Ricci and Riemann tensors corre- 
spondingly ||. 

In all physical applications the number of zero-modes does not change and we 
can drop the term with L(A + E). Returning to the initial variables and performing 
inverse Wick rotation, all curvature tensors and the metric are written in the "phys- 
ical" variables (i.e. without tildes). Moreover ln(64p Z 3 ) ^> ln(p) so that we can 
substitute ln(64pZ 3 ) instead of ln(64p / 3 ). Summarizing, we obtain an expression for 
the first quantum correction to the effective action: 



V x = -j dxdr v^ln(64p/ 3 )^ (A + E) , 



(6) 



with the metric 



1/2 f 1/2 



and 



IW = p L/2 f 



Vi = 2diip , / 



-p/f + v 2 
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Equation @ is the central equation of the paper. 

Miminizing the effective action (Bfi with the correction 
motion for the condensate can be derived as: 



^ - (Vy?) 2 - v + p - 21 p + ^ ln(64p/ 3 ) 



), the equations of 
dg^SVo , dESVo 



dp 5g^ + dp 5E 



, 
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8(p(y) 5g fiU (x) dr \dip 5E 



o 
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. 



The equations substitute for the Gross-Pitaevskii equation and generalize Eqns (|],[|) 
to account for the interaction of the Bose-condensate with the non-condensed frac- 
tion in the one-loop approximation. 

3. In conclusion, we have calculated the one- loop quantum correction for the 
effective action of the condensate. The correction comes from the interaction of the 
condensate with low-energy excited modes. To this end, we considered the excited 
modes as quantum particles moving in curved space whose metric (gravity field) 
was defined by the density and the velocity of the condensate. The equations of 
motion for the condensate are then analogous to the equations for the gravitational 
field where the right hand side (the energy-momentum tensor) is due to the presence 
of matter fields. The analysis does not depend on details of the external potential 
and can be used for variety of problems. The only simplification is the common 
hydrodynamic approximation. 
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